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A system and method for increasing 
jamming immunity in a GPS/INS system. 
The system includes a signal processor for 
receiving and down-converting GPS sig- 
nals to a baseband frequency to provide I 
and Q (sine and cosine) signals. An inte- 
grate and dump operation is performed on 
the I and Q signals to provide GPS range 
and range rate residuals which are uncor- 
rected from sample to sample. A Kalman 
filter estimates navigation state corrections 
from the GPS range and range rate resid- 
uals. A NAV function; responsive to the 
output of the Kalman filter and range er- 

ror and velocity errors signals, provides guidance information and signals for correcting at least one of earth centered earth fixed (ECEF) 
position, velocity, initial attitude and IMU alignment states which are dynamically calculated from the range and range rate measurements. 
A line of sight geometry function maps position and velocity vector information from the NAV function into geometric range and range 
rate scalar information and feeds the scalar information to the signal processor and circuitry to perform an integrate and dump operation. 
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INCREASED JAMMING IMMUNITY BY 
OPTIMIZING PROCESSING GAIN FOR GPS/INS SYSTEMS 

FIELD OF THE INVENTION 

This invention relates to a method and system for 
improving immunity to jamming in systems utilizing the 
global positioning system (GPS) . 

5 

BRIEF DESCRIPTION OF THE PRIOR ART 

Operations utilizing signals received from the GPS 
have been universally utilized both in conjunction with 
commercial and military applications. GPS and inertial 
10 measurement units (IMUs) have been combined or coupled 

together in the prior art to provide more effective 
navigation with the GPS correcting the IMU. This 
combination provides a synergistic result in that the 

effective bandwidth, of . the. .system .,.can be reduced, in an 

15 optimum fashion to provide improved^ t tracking. 

A problem with the use of the GPS is that it can 
easily be jammed. Systems have been devised in the prior 
" art to minimize the detrimental effects*: of GPS jamming. 
However, jamming and the possibility of ^jamming continue 
20 and still can ; present navigations - and ~ other : problems. 

: Accordingly , _\ systems ^ which provides.' increased immunity 
against such jamming are highly desirable and, in fact, 
, v . necessary in many . applications . c . ^It follows that any 
^improvement in "immunity against jamming by further 
25 V- reduction in bandwidth : is highly desirable. 

Traditional .~ GPS receivjeris^ 7;,!cbntain intermediate 
tracking loops which operate well at -high carrier to noise 
density (C/NO) (>23 dB-Hz) but: break down at low. C/NO 
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conditions.-' Delay lock loops (DLLs) are used for tracking 
code phase and phase lock loops (Costas loops) 'are used to 
track carrier phase. The outputs for- each loop type are 
used to generate pseudo range and delta (A) range 
measurements : respectively, which are used by a Kalman 
filter in formation of range and delta (A) range residuals 
(error signals- input to the Kalman -filter) . The noise 
bandwidths for these loops are usually on the order of 
0.025 and 1 Hz for the DLL and Costas loops respectively . 

The following list enumerates the problems associated 
with the traditional solutions to the jamming problem: 

*l-. : * Under jamming conditions (with low- C/NO) , narrow 
tracking loop bandwidths are usually employed/ resulting in 
temporally "correlated noise which is a suboptimal solution 
when viewed from the point of view of the Kalman filter. 
These narrow bandwidths also produce a correlation effect 
between measurements and process noise which tends to have 
destabilizing effects on system performance. 

2.' As jamming increases, the carrier tracking loop 
breaks down somewhere in the C/NO neighborhood of 18 dB^Hz . 
The fundamental reason for the loop break downs* is that the 
signal to : noise ratio (SNR) into the (Costas) loop, which 
is the product of carrier power/noise density (C/NO) and 
the coherent integration- time ( (C/NO) • Ti = SNR which is 
the signal to- noise ratio prior 'to the squaring operation) . 
When this happens, ' the informatioh loss ^through traditional 
non-linear loop error discriminants, the arctan function 
being an - example, ' becomes ' prohibitive, resulting in 
virtually no restoring force to any -loop perturbation. The 
only way to recover signal to noise ratio (SNR) is 1 to 
ihtegfate : longer and use narrow ibandwidths, resulting in 
correlated measurement problems -described "in * (1.) above. 
The loss of carrier loop measurements is key -since these 
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measurements provide. the- _ most* . accurate source of 
information necessary to maintain .the: .INS system alignment 
(i . e .,. system* errors small). ■ • z 

3. . Traditional schemes that, employ intermediate 
tracking loops are^ somewhat, sluggish to ; .reacquire (and 
• acquire) the GPS signals following.; a- jamout period. This 
. .is because, a certain amount of time rmust i be allocated for 
the tracking loops to pull, 7 in ■ to ; lock. - This time is 
somewhat large for narrow tracking - loop ; bandwidths which 
use long . time constants. -, t; .. L ^. y 

intermediate , tracking loops , employ- fixecl gains which 
do not "adapt" to. time varying signal conditions (e.g., 
../jamming .or loss, of data ,-.tp : be ; > . used- f or . data wipe off 
:(DWO) ) . • Data wipe off is -a. method jused to. improve signal to 
noise ratio.. Since I , and Q. data input to the baseband 
- .algorithm is bi-phase - shif ty; keyed ... (BPSK). ;4 modulated at a 
50 : Hz ..rate, , ; the signal polarity can. therefore change at a 
50 Hz rate .-(e.g. Q *=-.± A ,-sin.-. e ; ) .Clearly,, when this 
. happens, coherently adding, I or Q-rSignal samples over time 
can result in signal cancellation.^ DWO is a technique 
which uses ;r a priori estimates - of the- .50,. Hz. data stream to 
remove, this effect, thus allowing coherent integration and 
yits -benefits,. DWO is . normally ^accomplished by ^demodulation 
of./ the data^ prior to . j ampul: - .and then . storing of, the data 
. ; f or later.- use . This; demodulation .must take place for the 
primary purpose of navigation/ . SQ : ,the ,availabili'ty of most 
of- the required data .is, a given.. , When, this happens, the 
tracking, loops, are in. >a .-state iOf .random walk with a 
..: .r di/yergence rate - determined by v the. ; f ixed gains;, employed. 

.This creates measurement,- editing and recovery problems 
• : • when the;- signal conditions, again -become favorable. This 
/ creates a-yery ^suboptimal situation, , in which data used for 
.-» y . p^O.is misestimated.. 
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The standard performance metrics when discussing GPS 
operation under jamming conditions are receiver threshold 
jamming to signal levels (J/S) and carrier to noise density 
(C/NO) . Under conditions of wideband jamming, C/NO (dB-Hz ) 
= 73 - J/S (dB) . Therefore, under conditions of high J/S, 
the signal to noise ratio (SNR) decreases proportionally 
for signals delivered to baseband signal processing 
algorithms' (intermediate baseband tracking loops). 

Ah approach to reducing tracking loop noise bandwidth 
in order to navigate during modest- amounts of jamming .is 
inertial measurement unit (IMU) -inertial aiding techniques 
of code and' carrier tracking loops* Assuming that a 
nulling antenna is not used, tracking loops currently phase 
lock and code lock to J/S levels of approximately 55 dB and 
60 dB respectively. The outputs of these loops produce 
range and delta (A) range measurements to a Kalman filter 
(KF) , generally an 18 state navigation Kalman filter though 
it should be understood that the number of states of the 
Kalman filter is optional. Above these jamming levels, 
complete loss of GPS measurements occurs, resulting in the 
inertial navigation system drifting in a divergent fashion. 

SUMMARY OF THE INVENTION 

In accordance with the present invention, a procedure 
is provided which uses integrate and dump techniques 
operating on I and Q data to directly produce residuals 
which are input to a Kalman filter used to correct 
navigation errors without the use of intermediate tracking 
loops. The techniques described provide optimal processing 
gain while retaining the Kalman filter * optimality 
requirements of uncorrelated measurement errors. This 
technique provides both code and carrier measurement at 
jamming levels heretofore unheard of . ; ■■ Qualify carrier 
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meas.urement at J/S- levels of 70 and 75 .dB are attained for 
. ■. carrier phase and frequency measurements respectively, 
given a joint standoff weapon ; (JSQW) quality IMU, which 
.uses a low grade (3 sigma) ring laser gyro. Code tracking 
5 through .the.. Kalman filter . .is also attained in ,the 

, neighborhood of 7 5 dB - J/S. . . The high .quality carrier 
.measurement obtained at these extreme . j $mming levels is an 
important feature of the present, invention. The accuracy 
of -the measurement .maintains . the . inertial, navigation system 
10 (INS) -alignment by : keeping-. navigation »■ errors small, 

particularly velocity, and .attitude., errors, , Furthermore, 
; i these techniques produce jamming, immunity without the use 
of , but not precluding, a ; nulling antenna. t This invention 
.therefore offers a solution to the .jamming problem. In 
;15.. -summary, in accordance ; with the present .invention, carrier 
and code tracking is possible . (assuming. ,DWO) t with 15 dB 
.- ; more jamming than the current state : of .the art. Without 
the. use of- DWO, a 6 dB improvement is, ^tq. be expected. „ 
. ■ Briefly, appropriate,, coherent, .and.- non-coherent 
20 . integration.- of I and ,Q data y-to. r minimize, nonlinear 
processing, losses which Is optionally combined . with bi- 
phase shift keyed (BPSK) data cancellation techniques known 
as data wipe off (DWO) produces* residuals . which, directly 
drive the .navigation .filter. . ..The . following discussion 
25--.- . describes-:. :,the. processing, for .-, three ; different .types- of 
. ;.; - ; i. •. residuals, all of which input I and Q- data ; at .a relatively 

fast. rate. . . - - 1 \ .:■ ■■ - \- .- u 

•< , -'-.GPS .signals are -received .and down-converted to a 

baseband frequency . and are .input,. to a -signal processing 
30 ,- chip , (SPC) ^.in standard . manner . ... The SPC. . contains a 
plurality of GPS channels to^ enable . tracking of a plurality 
- . ...-of , satellites simultaneously. ..Accordingly, signals from 
: * : - plural^ satellites ; are received and processed in the SPC in 
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parallel to deliver standard samples of I and Q signals 
(sines and cosines) derived from the GPS system signals 
down 'converted to the baseband software. 

In the prior art, intermediate tracking loops 
utilizing 1 & Q signals received from the SPC were used to 
track range and A range (range rate) to each satellite in 
order to calculate a residual (error signal that drives the 
Kalman Filter) . This residual is sent to a Kalman filter 
(KF) which then sends correction data to a generic 
navigation (NAV) function. A strap down NAV function is 
preferred, but is optional. Based upon the signals 
received from the IMU and the KF, the navigation function 
corrects earth centered earth fixed (ECEF) ; position and 
velocity as well as initial attitude and < IMU alignment 
states which are dynamically calculated from range and 
range rate measurements. The NAV function is standard 
known software^ capable ; of blending raw IMU data with 
navigation corrections derived from GPS measurements to 
produce a corrected navigation solution at very fast update 
rates. These NAV update rates are generally much faster 
than the GPS measurement rate, but the only requirement is 
that the updates "occur at least as fast as the measurement 
rate of occurrence. The NAV function then feeds back an 
estimate of the range - and range rate -to each satellite 
through a line of sight (LOS) geometry function to the 
tracking loops which then tracks' out the aiding errors. 
The' LOS geometry function maps navigation position and 
velocity : (vector) information into geometric range and 
range rate - (scalar ) information for a given user and each 
satellite tracked. • These estimates are input to the SPC 
and I & D blocks, based on the system's navigation 
solution, as- range and range rate 'for' each' " satellite 
tracked. These estimates are used in a : standard' manner to 



^WQ 98/48295 PCTAJS98/04654 

7 



produce (or. measure) -the input ^1 & Q values which are 
input, for example, at a. 200 Hz. rate. The aiding errors 
are the errors in the said range and range, rate estimates 
described above . , The Kalman filter, by. making corrections 
5. . to all of the navigation : states ; (position, velocity, 
. -attitude, etc.)., consequently also reduces the aiding 
, errors to zero (i.e.., producing resi duals , which are on the 
. average zero) , thus stabilizing. ; the system. . 

By ;the present invention, the intermediate tracking 
10 . ; loops are- replaced, by an; .integrate -and. dump (I & D) 
operation which places, the- residuals ciirectly into the 
Kalman; filter. - By,- .this method, :::? each r measurement is 
independent from 1 sample to , sample - and .loop closure is 
..-accomplished by : the Kalman filter v (and navigation (NAV) 
5 „ ^function) > as opposed to- ,the tracking^ loops which are 
themselves , closed . loop systems^, r _ This .procedure takes 
advantage of certain mathematical/statistical properties of 
the residual estimates used, to, • enhance- SNR and thus 
provides a more stable, system at -higher jamming levels. 
2.0 The techniques described provide, -optimal processing 

gain, while retaining., the. ... .Kalmari; • filter optimality 
.requirements of uncorrelated -measurement errors. Each 
integrate and dump (I : &...D-) algorittinu -inputs the I and Q 
. signals from- the SPC, which v arp ninety degrees out of phase 
25. with.,, each other, ,and ; received , at a, predetermined rate. 

These signals, are .added, up -coherently (.in. phase) over a 
designated,- time frame (Ti) v . This : , time- -frame- is determined 
as -a function of LOS- range- rate .covariance, calculated from 
Kalman filter covariances in order ; -tp determine the optimum 
30 predetection integration ...time T x . ■- This functional (T-i is 
: ... v>- approximately .1 /range rate uncertainty) relationship can be 
r , : . implemented as an equation qr. ,as a - look-up table. The 
i v outputs-' .from 7 each coherent- summation . £ ( ) , cross and dot 
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products are formed at the Ti rate and the cross and dot 
products are then summed over a different time interval of 
length Tk. This -time, Tk, which is chosen to a targeted or 
desired signal to noise ratio, is determined by SNR 
estimate's computed by- the system to provide a pair of 
second summations or integrals, one for the dot product 
signal and- one for the cross product signal. Then an 
arctangent (also set forth herein as arctan or atan) 
function- of the summed cross products and dot products 
[atan2 (cross, dot)] is taken and divided by the integration 
time (T-i) to- obtain the range rate residual directly. This 
residual' f 'is then input to a Kalman filter (KF) 

The algorithm which computes range residuals, receives 
the I and Q signals from the SPC at a predetermined rate. 
These signals- are added up coherently over a designated 
time frame (Ti) for each of the I and Q signals 
individually to provide summations or integrals for each of 
the I and Q signals. Such procedure is provided for both 
an early channel and a' 1 late channel. In each channel, the 
■I signals and Q signals are squared, this being an envelope 
detect type of operation". The squared signals are summed 
-for the time' period Tk. 

The algorithm which computes a phase lock loop (PLL) 
range residual is a biased estimate in contrast to the 
range/ residual' described above which provides an unbiased 
estimate of the range to the satellite but is quite noisy. 
The Kalman filter separates the noise from the signal 
optimally as long as it is advised of the amount of noise 
present in the ;total- signal sent thereto. The more noise 
present, the slower the filter operates. Here, the Kalman 
filter is' used * to close a- phase lock .loop (PLL) by 
maintaining all of the energy in the I (cosine) channel and 
providing no . energy in the Q (sine) . channel -or - keeping that 
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signal at zero. Although the PLL .provides a -measurement 
v . with; two orders of magnitude less noise, a problem with the 
PLL is. that 2n looks like 4n . which., looks .like 6n, etc. and 
has -ambiguity. Each 2n represents . -a 0.2. meter bias. 
5 : . Accordingly, . the . system : is allowed .rto converge using the 
unbiased estimate and then,, when-. £he ; error is low, the PLL 
mode is used. In accordance with this a-lgori-thm, the I and 
Q signals from the SPG- are -.received at ,.a. predetermined 
rate. ; These signals. , are, ,-added ;.;up^ v coherently over a 
0- .designated time,, frame (Ti) for each . of ? -the I- and Q signals 
individually to provide. summations^ or inte.grvals .for each of 
the I and Q signals The arctangents :; -of the , s t ummed signals 
are ; - taken to provide- a range residual to, the;'Kalman filter. 
; ^Although this estimate ■ provides a ; modulo^ 2n bias in range, 
5 : . it provides the most accurate measurement, with regard to 
the estimation of velocity -and attdtude , when compared to 
.the other residual described, herein v 

BRIEF DESCRIPTION OF THE DRAWINGS ^ : 

0 ■* , FIGURE - 1 is- a -block-diagram of -/an-architecture for use 

. •in conjunction with the pr.esentv invention ; - : .-.*, 

FIGURE 2 is a block diagram, representing the algorithm 
, for. the integrate . and dump:. (I ,&.-*. D>* block of FIGURE 1 for 
range -rate residual algorithms;^.. '/■-•:.;*> ; • - 
5 . FIGURE 3 is. a block diagram; representing the algorithm 

v for .early, late * range ; residuals -which /replaces : the I & D 
block .-of FIGURE 1. ; and-;.,:-;;... 

r. FIGURE 4 is*. a block. diagram, representing the- algorithm 
• v . forr-the .phase lock loop -(PLL).' /range residual-. • .* 

DESCRIPTION OF THE PREFERRED EMBODIMENT . 

■ rl: ; ^Referring .to FIGURE:, 1, .thefe. is; shown- ;a bloqk diagram 
|t ;:. ; ,:'of-'.anc -architecture .- f or - uset ini conjunction with, .the present 
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invention. GPS signals are received and down-converted to 
a baseband frequency by a standard signal processing chip 
(SPC) in standard manner. The SPC contains a plurality of 
GPS channels to enable tracking of a plurality of 
satellites - simultaneously. Accordingly, signals from 
plural satellites are received and processed in the SPC in 
parallel to deliver standard I and Q sampled signals (sines 
and cosines) . derived .from the GPS system signals to the 
baseband, software;- v 

* The intermediate tracking loops of the. prior art are 
replaced by. an integrate and dump (I & D) operation which 
inputs the residuals directly into the KF. 

The- algorithm, uses integrate and dump (I & D) 
techniques operating on the I & Q data to directly produce 
residuals input to, for, example,, an 18 state KF without the 
use of ' intermediate ■. tracking loops. The technique 

■ described provides optimal processing gain while retaining 
the KF optimality requirements of uncorrelated measurement 
errors. This technique provides both code and carrier 
measurements .at jamming levels heretofore not. attainable . 
Quality carrier loop measurement at J/S levels of 70 and 75 
are attained , for carrier phase and frequency measurements 
respectively, given a JSOW quality IMU.: The high quality 
carrier measurements at these extreme, jamming . levels are an 
important feature. of the invention. The accuracy of these 
measurements, maintains,, the INS . - system errors small, 
particularly velocity .and. attitude errors. ' Furthermore, 
these techniques produce j amming : immunity • without the use 
of, but not. precluding,, a nulling, antenna, 

Referring to FIGURE .2, .there is shown a block diagram 
representing the algorithm, which is contained in the 
integrate and . dump ( I & D) block of FIGURE 1- for the range 
rate residual algorithm. The I and Q signals from the SPC, 
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which ninety degrees out of phase' with". each* other, are 
received : at a predetermined rate, this rate being shown as 
; - - 200 Hz ' in the preferred embodiment v , * These signals are 
■ added up coherently (in phase) over a designated time frame 
5 ■ (Ti) for each of the .< I and- Q ' signals individually to 
provide* a summation, (or integral: (£( )0 )' for each of the I 
and Q signals over the time period:: Ti . From these 
> summations- £() , cross and' dot ^products- are* formed at the Ti 
rate and the cross and dot products ^ar.e, then summed over a 
10 different -time interval of length Tk to provide a pair, of 

' second summations or integrals, - one; for ifehe. cross product 
signal, and one for. the /.dot : . product ..signal . s Then an 
arctangent function of the* summed-vcross ^products and dot 

products- is taken and divided -by: the-: integration time (Ti) 

15 - to obtain- the ■ residual -directly . . This: jresidual is then 
.vv - input to the Kalman filter (KF) Van FIGURE : 1. 

Referring to FIGURE -3,.. there - isv shown a block diagram 
"representing the algorithm' for early elate range residuals 
which represent the P & D . block of iiFIGURE . 1 . . This range 
-20 - residual .provides a globally/ stable: unbiased estimate of 
, v the range^ to the satellite ^but is quite noisy (i^e., large 
standard deviation) .'• In. accordance- with this . algorithm, 
the I and Q -signals from the: ^SPG . rare, - received at a 
- . predetermined rate, this -rate obeing vshown as*. 200 Hz in the 

•:25 preferred" V embodiment.. These.- -signals are;. added ! up 

: coherently -over a designated -time frame /;.(Ti) for each of 
the I and Qf signals individually . to ^provide, .summations or 
integrals for .each of; the :;I- and; Q.\signals .As .can be seen 
from FIGURE ,3.,. such, procedure Ls 'provided, for- each of early 
'.30 channel and late channel '(/separated by ± one half chip) . 

: .-In each of, early:- channel- and.- late channel,. * the I signals 

})s r.i and Q signals are squared, this, being an envelope detect 
i:;. r . type ;of, operation . _ The squared :signals are sjummed for the 
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time period Ti. Envelope detection then takes place for 
each channel prior to computation of the range residual 
which is the : difference between early channel and late 
channel divided by the sum of early channel and late 
5 channel, an equation therefore being shown herein below. 

Referring- now to FIGURE 4, there is shown a block 
diagram representing, the algorithm for the phase lock loop 
■ (PLL) range residual. The Kalman filter can separate the 
noise from the signal optimally, as long as it is advised of 

10 the amount of noise present in the residual sent to the 

Kalman filter. The more noise present, the slower the 
filter operates. Here, a phase lock loop (PLL) is used to 
maintain all of the energy in the I channel and provide no 
energy in ,the Q (sine) channel or keep that signal at zero. 

15 A problem with the PLL is that 2n looks like 4n which looks 

like 6n, etc. to provide an ambiguity. This residual 
provides local stability at increments of 2n, but each 2n 
represents a 0.2 meter bias. Accordingly, the system is 
allowed to converge using the unbiased estimate (as 
; 20 described with reference to FIGURE 2) and then, when the 

. system converges, the mode of FIGURE 4 is used. 
Convergence is determined by checking the Kalman filter 
covariance. - In accordance with this algorithm, the I and 
Q signals from the SPC are received at a predetermined 

25 rate, this rate being shown as 200 Hz in the preferred 

embodiment. These signals are added up coherently over a 
designated time frame .(Ti) for each of the I and Q signals 
individually to provide summations or integrals for each of 
the I and Q signals. The arctangents of .the summed signals 

30 are taken to provide the output which is the input to the 

Kalman filter as a range residual. 

The following is a description of an algorithm used 
herein in conjunction with the present invention which is 
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to be used in • conjunction with the drawings as described 
hereinabove. 

The range and range .rate aiding is utilized to form I 
and. Q data which is mathematically : modeled as the sines and 
5 cosines of .the -carrier phase error as follows: . 

• I *J= D(t:);*: DWO estim ate (t) • R XX (t,) COS (B e ( U ) . 

: Q = D(t)- DWO estimate (t) • R X x(r c )sin(© e (t) ) • 
" where D*(t) is the 5G Hz rate, R xx lis -the PN autocorrelation 
• triangle 1 ; x; is the range- error in chips (1 chip = 30 
0 meters) and B e is the range'error in radians. The carrier 

phase error B e (t) represents the range : error in radians 
: referenced to the carrier frequency (I. e.' ; , LI = 1.54 .GHz 
or ; 154f 0/ : where fO = 10.23 MHz ) This- range error is the 
error between the true range to a given satellite and the 
'5 range estimated from navigation state /variables (these 

: include ECEF position and velocity as well as' user time 
bias and bias rate variables ). ' 

The" code phase mis-tuning T c , in ■ cycles @ fO, also 
represents the iange error and' has the" effect of a 
0 ' correlation loss R xx on the ^signal. Note that R xx denotes 
the standard auto-correlation triangle for a PN code. In 
the case of the P-code usecT for GPS , this triangle spans 
(±) one p"chip J (cycle @ f Q or about 30 meters) . Not shown 
aire frequency mistuning fosses " (sinx/x) and additive 
5 - Gaussian noise for each I and Q 'data pair'. Note also that 
D(t) rtiodulates the carrier at a 50 Hz rate '"and' is removed 
(or wiped" off) 1 with estimates of D(t) which are passed to 
this routine : from a data demodulation function'. Most of 
this modulated data is required to ! be known independent of 
0 ; signal processing benefits, 'since the data' contains 
ephemeris information required for navigation purposes. 
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In order to estimate pseudo range and range rate 
residuals input to the navigation Kalman filter, the 
following calculations are made: 

For the first residual, first the I and Q values are 
5 coherently integrated for a specified period of time which 

depends upon the range rate uncertainty (LOS covariance) 
calculated by the Kalman filter according to: 

N 
i = \ 

N 



The bars over the . I and Q data indicate that they are 
average orthogonal : components of a phasor . The argument 

10 - ~ for each component represents, with good approximation, the 
average phase error evaluated at . the midpoint of each 
coherent integration dwell. 

Using the above coherently derived I and Q values, 
three different residuals are formed to drive the eighteen 

15 state extended -Kalman. filter. FIGURE 1 illustrates the 

high level structure for each of these residuals, which are 
produced as integrate and dump (I & D) residuals to the 
Kalman filter . The*- three residual types are described as 
follows: .<■'.. *- 

20 By coherently integrating I . and Q values for both 

' early and late- -P-^code state : setups, the, .unbiased pseudo- 
range residual - is -formed, as. follows from coherent samples 
according to FIGURE 2 : \ - ■ . 
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E - [(t ai E ) 2 * (Q E ) 2 )ft m : 

l = [(tuT L ) 2 '+ (Q]) 2 )i)V n ■ ■-" 

The range residual is then formed as: 

E+L 

- G is a gain factor derived from C/NO estimates to 
compensate for gain depression, effects . This gain factor 
is contained in a lookup table and is a function of C/NO. 
5 ; : . The residual 5R, scaled to meters ,.. .then drives the 

extended Kalman filter every- t k - seconds ... -This residual is 
an- unbiased range error., measurement ..from which the 
navigation ECEF position., and: clock; .states are directly 
observable. This residual: makes, .use of global properties 
10 o'f ■. p-code that produces- zero-prompt channel correlation 

— loss .when and only when) the pseudo-range residual is zero 
(R xx (0)+1) or, alternatively^.when E=I*.. ... rUsing this present 
. invention, the Kalman . filter, has.. - taken . the place of 
. . traditional tracking loop filters in; order to close a delay 
15 . lock loop.; i The taction of -this closed, loop system is to 
estimate navigation states that directly drive the pseudo- 
range residual .to zero.. - ... > / 

' -. The second .residual, estimated .to -directly drive the 
Kalman filter is an unbiased , estimate^ of pseudo range rate. 
20 This residual has as inputs coherently; generated I and Q 

data sampled from the prompt code correlator. The pseudo 
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range rate residual is calculated as follows as shown in 
FIGURE 3: 

SR = dWjdt = ATanlicross, dot)ldt 

i = i 

'' = i 



where I ' and Q data are coherently averaged over adjacent 
time steps (h and n-1), separated in time by , dt . This 

5 residual : 'is then scaled from units of radians @ LI to m/s 

and then directly input to the Kalman filter. What has 
happened in this case is that the Kalman filter has been 
used in : the plate of traditional tracking loop filters in 
order to close' l a frequency lock loop. This measurement 

0 directly couples to- the velocity and time bias rate state 

variables of the Kalman filter. 

The ' third - residual for direct .: use into the Kalman 
filter is treated as : a range residual and : provides a 
-locally stable but biased range residual (by -modulo 2n or 

5 0.195 meters @ LI). This residual is calculated as follows 

as shown in- -'FIGURE- 4: 

(a.) Si? = ® e (t mid ) = £ (4RCTAN2(QJ)) i 
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(b.) 6R = ® e (t mid ) = Y. (ARCTAN(Q/1)X± . 

1=1 



In case (a..), the arctan2 function is used when DWO is in 
use. When DWO is not in use) then the two quadrant (arctan 
(y,x)) in case (b.) supra is used. The choice is 
application specific. 
5 The third residual /represents the pseudo-range error 

(modulo 2n) with units of radians @ the LI carrier 
frequency. In contrast, the first residual described above 
is an unbiased estimate of pseudo range to each satellite 
> but @ fO. -When both residuals are scaled to ^meters, the 
10 r biased estimate in the third residual : is ;j inherently 154 
. -times. more accurate than the, first residual . _ This inherent 
accuracy difference produces.; carrier^ measurements into the 
Kalman filter with measurement noise on the order of 10 to 
' .< 20 mm (one sigma). at C/NO values ; in^ the - neighborhood of 
15 , 2*. dB-Hz. These measurements, .produce. ; the most accurate INS 
alignment . at - low. C/NO , values; t.apgroaching . 2 dB-Hz) by 
coupling to high derivative /.states, , through the Kalman 
... filter dynamics equations . . What has; happened in this case 
- : - is -that ,;.the- Kalman. filter, has; :. been : used in place of 
20 .', traditional! tracking loop fil-ters A in order to close a phase 

■ ' - r . lock- loop. ... I : ■ . : ; -'. • " " 

In order to minimize the. navigation position bias for 
this estimate, the system first converges using the first 
two ^.residuals; : Using. - Kalman filter error covariance 
25 estimates as indicators of system convergence, the third 

range residual is then used for ultimate INS alignment 
accuracy at very low C/NO values. The second range rate 
residual may be used simultaneously with the first and 
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third residuals. It should.be noted that/ for all three 
residuals, the Kalman filter is iterated every K*N*dt 
seconds. Due to throughput limitations, non-coherent 
integration in the third- residual may be particularly 

5 useful* when applying this method at high C/NO conditions 

(option (b.) above). 

Though the invention has been described with reference 
to a specific preferred embodiment thereof, many variations 
and modifications will immediately become apparent to those 

0 • skilled in the art. -It is therefore the intention that the 
appended claims be interpreted as broadly as possible in 
view of "the prior art to include all such variations and 
modifications. * 
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CLAIMS - . ~ 

1 . .. A system for increasing j amming . immunity in a 
• GPS/ INS system which .comprises : 

(a) a signal processor ( for receiving and down- 
5 ' converting. 'GPS signals to a baseband -frequency to provide 

I and Q (sine and cosine) signals; ,, .■ . 

(b) means to- perform an. Integrate and dump operation 
- - on. said I and Q signals to provide -;GPS. range and- range rate 

residuals which are uncorrelated from sample to sample; 
0 ; (c) a Kalman filter for estimating navigation state 

. corrections from said GPS range and- ranges-rate, residuals; 

. (d) . NAV means responsive to the output of- said Kalman 
-filter and range error and velocity ...error; signals to 
provide guidance information and signals for correcting at 
5 least one of earth centered earth fixed (ECEF) position, 

velocity, initial attitude and IMU alignment states which 
are dynamically calculated from range and range rate 
measurements; and 

(e) line of sight means for mapping position and 
0 velocity vector information from said NAV means into 

geometric range and range rate scalar information and 
feeding said scalar information to said signal processor 
and said means to perform an integrate and dump operation. 

5 2. The system of claim 1 wherein said signal 

processor includes a plurality of GPS channels to enable 
tracking of a plurality of satellites simultaneously. 
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3. 'The system of claim 1 wherein said means to 
integrate and dump includes coherently, adding the I signals 
and the Q signals separately over a predetermined time 
period (Ti )' to ; provide a summation for each of the I and Q 
signals over the time period Ti, forming from these 
summations cross and dot products at the Ti rate, summing 
the cross* and dot products over a different time interval 
of length Tk to provide a pair of second summations, one 
for the cross product signal and one for the. dot product 
signal, and taking an arctangent function of the summed 
cross products and dot products and dividing the arctangent 
function by the integration time (Ti) to obtain said 
residual directly. - 

4. The system of claim 2 wherein said means to 
integrate and dump includes coherently adding the, I signals 
and the Q signals separately over a predetermined time 
period (Ti) to provide a summation for each of the I and Q 
signals over the time period Ti, forming from these 
summations cross and dot- products at the Ti rate, summing 
the cross and dot products over a different time interval 
of length Tk to provide a pair of second summations, one 
for the cross product- signal and one for the dot product 
signal and taking an arctangent function of the summed 
cross products and dot products and- dividing the arctangent 
function by the integration -time (Ti) : to obtain said 
residual directly. - " 
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5. The system of claim 1 wherein said means to 
integrate and dump includes an adder .for coherently adding 
the I signals and the Q signals separately over a 
predetermined time period (Ti ) :to ; provide a. summation for 
each of -the I and . Q signals, over the time period Ti for 
each of early channel and late- channel, , squaring circuitry 
for squaring the received I and Q ..signals -in- each of early 
channel >and" late channel to provide an envelope detect, a 
summer for summing the squared signals separately in each 

0 of the early channel and late channel .for .the time period 

\.Ti, . an ■ envelope- detection circuit providing-, envelope 
, detection for each channel -.-prior;..; to -computation of ^ the 
range residual, said envelope detection being circuitry to 
provide the difference between early channel and late 
5- ■- channel .divided by the. sum of ...early channel and late 
channel. . . ; • 

6. The system of claim: 2 wherein ; said means to 
^integrate: and dump . includes . an adder., for coherently adding 

0 : . .the I- signals and the - Q,. signals, separately over a 
predetermined time-period ; (Ti) to_, provide a- summation for 
.each ;of;the I and Q ■ signals, over- the time period Ti for 
each, of early channel and late-channel, squaring circuitry 
for- squaring.- the received. I- and t Q -signals in each of early 
5 channel andilate channel-to. provide ^an ; envelope . detect, a 

.summer for summing the^ squared signals separately in each 
of the early channel and late channel,- for the time period 
Ti, an envelope detection circuit providing envelope 
detection for each channel prior to computation of the 
0 range residual, said envelope detection being circuitry to 

provide the difference between early channel and late 
channel divided by the sum of early channel and late 
channel . 



WO 98/48295 



22 



PCT/IJS98/04654 



7. ' The system of claim 3, further including, 
responsive to said range residual reaching a predetermined 
value, substituting for said means to integrate and dump 
means for receiving the I and Q signals from the SPC at a 

5 predetermined rate, an adder for separately coherently 

adding said' I and' Q signals over a designated time frame 
(Ti) for each of the I and Q signals individually to 
provide summations for each of the I and Q signals and 
means to take the arctangents of the summed signals to 
10 provide" the output which is the input to the Kalman filter 

as a range residual. 

8. The system of claim 4, further including, 
responsive to said range residual reaching a predetermined 

15 value, substituting for said means to integrate and dump 

means for receiving the I and Q signals from the SPC at a 
predetermined rate, an adder for separately coherently 
adding said I and Q signals over a designated time frame 
(Ti) for each of the I and Q signals individually to 

20 "provide summations for each of the I and Q signals and 
means to take the arctangents of the summed signals to 
provide the output which is the input to the Kalman filter 
as a range residual. - " ■ : 
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; - , 9. The system of claim 5, further including, 

responsive to said range residual .reaching a predetermined 
, value, substituting for said means, .to integrate and dump 

means for receiving ..the . I and - Q . signals . from _ the. SPC at a 
5 ; predetermined rate, - an adder . for separately coherently 
adding said I and Q signals over, a -designated time frame 
(Ti) for each of the I and Q . signals. ( individually to 
provide summations for each of . the f I^and Q signals and 
means to take the • arctangents of the summed signals to 
10 • provide the output which is. the input to, the Kalman filter 
as a range residual. f . ■ - ; . : - .*": -, 

... ; 10., The system of claim 6,.,. further including, 
responsive to said range residual .reaching a . predetermined 
, - 1-5 . value, substituting for said .means to : . integrate and dump 
means for receiving /the I and . Q, signals, from tl^e SPC at a 
predetermined rate, an adder, for . : separately coherently 
adding said I and Q signals over . a ~ designated time frame 
; (Ti) for each of the I and Q. signals , individually to 
• 20 provide .summations for each of. .the if ,I ..and Q signals and 

means ; , to take the arctangents ;i qf ..the summed signals to 
provide the output which is the ..input to the Kalman filter 
as a range residual. 
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11'. A method for increasing jamming immunity in a 
GPS/INS system which comprises the steps of: 

(a) receiving and down-converting GPS signals to a 
baseband frequency to provide I and Q (sine and cosine) 
signals; 

(b) performing an- integrate and dump operation on 
said I and Q signals to provide residuals which are 
uncorrelated from* sample to sample; 

(c) performing * a Kalman filtering function for 
estimating navigation .state corrections from said GPS range 
and range rate residuals; 

(d) responsive to the output of said Kalman filter 
and range error and velocity error signals , providing a NAV 
function to provide guidance information and signals for 
correcting at least .one of earth centered earth fixed 
(ECEF) position, velocity, -initial attitude and IMU 
alignment states which are dynamically calculated from 
range and range rate measurements; and 

(e) mapping position and velocity vector information 
from said NAV -function into geometric range and range rate 
scalar information and feeding, said scalar information to 
said signal processor and said < means to. perform an 
integrate and dump operation. 

12. The method of claim 11 .wherein said signal 
processor is provided with a plurality of GPS channels to 
enable tracking of a plurality of satellites 
simultaneously . 
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- 13. ' The method of claim 11. wherein said integrate and 
dump operation includes.. coherently adding the, I signals and 
the Q signals separately over a • predetermined time period 
(Ti) to provide a ■summation for each* of /the r and Q signals 

5 over the time period Ti, forming from these summations 

cross and dot products; at- the , >Ti, rate,: summing the cross 
and dot products over a different time ^interval of length 
Tk to provide a pair of , second , summations, / one for the 
• cross product signal and one for the dot, product signal and 

0 - taking an- arctangent function- of 'the summed: cross products 
and dot products and dividing ..the arctangent function by 

... * . the integration time. (Ti) to obtains said residual directly* 

14. . The: method of claim 12 wherein - said Integrate and 
5 L dump operation includes coherently adding^ t lie .1 signals and 
V* the Q : signals separately. -.over a predetermined time period 
(Ti) to provide a summation 1 for , each of the I and Q signals 
over the time period Ti,.-, .forming ,.from these summations 
cross and dot products at the Ti. rate, summing the cross 
0- and dot products over a: dif ferent; time/, interval .of length 

Tk to provide a pair of s.econdj ^summations, one for the 
- -i cross product, signal and one .for; the dot . product signal and 
taking an arctangent functionrof the .summed cross products 
and dot products and dividing the arctangent function by 
5 '- the integration time-. (Ti) Ho obtain: said residual directly. 
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15/ The method of claim 11 wherein said integrate and 
dump operation includes coherently adding the I -signals and 
the Q signals separately over a predetermined time period 
(Ti) to provide a summation for each of the I and Q signals 
5 over the time period Ti for each' of early channel and late 

channel, squaring the received I and Q signals in each of 
early channel and late channel to provide an envelope 
detect, summing the squared signals separately in each of 
the early channel and late channel for the time period Ti, 
10 providing envelope detection for each channel prior to 

computation of the range residual, said envelope detection 
being the difference between early channel and late channel 
divided by the sum of early channel and late channel . 

16. The method of claim 12 wherein said integrate and 
dump operation includes coherently adding the I signals and 
the Q signals separately 1 over a -predetermined time period 
(Ti) to provide a summation for each of the I and Q signals 
over the time period Ti for each of early channel and late 
channel, " squaring the received I and-Q" signals in each of 
early channel and late channel to provide an envelope 
detect, summing the squared signals separately in each of 
the early channel and late channel for the time period Ti, 
providing envelope ' detection for ' each channel prior to 
computation of the range residual, said envelope detection 
being the difference between early channel and late channel 
divided by the sum of - early channel and late channel. 
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17. The . method ... of-, claim -13, further including, 
responsive to said range residual . reaching a predetermined 
value, substituting for said .integrate and dump operation 
receiving the : I. and Q signals vfrom the SPC at a 
5. predetermined, rate, separately coherently^ adding said I and 

Q signals over a designated time frame .(Ti) for each of the 
I and.Q signals individually - to provide summations for each 
of the I : and Q signals and- taking, the ..arctangents of the 
.summed, signals to -provide the output: which; is the input to 
0 the Kalman filter -as a range residual . s ;; ; - , 

18^ The method .of . claim 14> t further , including, 
-responsive to said range rjssidual reaching a -predetermined 
value, substituting for said integrate and dump operation 
5 • .receiving the I and . Q .signals., .from the SPC at a 
. . predetermined rate, separately^ coherently, adding said I and 
Q signals over a designated time frame . (Ti) - for each . .of the 
I and \.Q signals individually to. -provide summations for each 
of the I and Q . signals ( and -taking the. : arctangents,<.of the 
0 - summed signals to provide „ the output .which ; is the input to 
. . the Kalman filter as ,a range residual.-.. 

,...19. .-, The , method of . claim 15, .; further including, 
responsive to - said range residual. reaching a predetermined 
5. value, substituting- f or . said integrate : and dump operation 
... .receiving ,.the I - and .Q- .signals vfrom the SPC at a 
, : , predetermined rate, separately, .coherently adding said I and 
Q signals over a designated time frame (Ti) for each of the 
I and Q signals individually to provide summations for each 
0 of the I and Q signals and taking the arctangents of the 

siammed signals to provide the output which is the input to 
the Kalman filter as a range residual. 
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20. The method of claim 16, further including, 
responsive to said range residual reaching a predetermined 
value, substituting for said integrate and dump operation 
receiving the I and Q signals from the SPC at a 
5 predetermined rate, separately coherently adding said I and 

Q signals over a designated time frame (Ti) for each of the 
I and Q signals individually to provide summations for each 
of the I and Q signals and taking the arctangents of the 
summed signals to- provide the output which is the input to 
10 the Kalman filter as a range residual. 
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